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We study a dynamic process of disentanglement by considering the time evolution of bound 
entanglement for a quantum open system, two qutrits coupling to a common environment. Here, 
the initial quantum correlations of the two qutrits are characterized by the bound entanglement. In 
order to show the universality of the role of environment on bound entanglement, both bosonic and 
spin environments are considered. We found that the bound entanglement displays collapses and 
revivals, and it can be stable against small temperature and time change. The thermal fluctuation 
effects on bound entanglement are also considered. 

PACS numbers: 05.30.-d, 03.65.Ud, 75.10.Jm 



I. INTRODUCTION 

Entanglement [1], as an essential feature of quantum 
mechanics, helps us to distinguish the classical and quan- 
tum nature of matter world. It has become a key ingredi- 
ent in quantum information processing, such as quantum 
computing, quantum teleportation and quantum cryp- 
tography [2]- [4]. On the other hand, generally a realis- 
tic system is surrounded by an environment. Thus the 
effects of the quantum decoherence such as quantum de- 
phasing on quantum entanglement should be considered 
for quantum opene systems. It is reasonable that when 
we study quantum effects induced by entanglement, the 
two-particle system should hold phase relations between 
the components of the entangled states. Thus perceiv- 
ably, due to the interactions with environment, we can 
expect that the dephasing of two-particle system can 
demonstrate some exotic properties. 

Recently, Yu and Eberly [5] showed that two entan- 
gled qubits became completely disentangled in a finite 
time under the influence of pure vacuum noise. Surpris- 
ingly, they found that the behaviors of local decoher- 
ence is different from the spontaneous disentanglement. 
The decoherence effects take an infinite time evolution 
under the influence of vacuum noise while the entangle- 
ment displays a "sudden death" in a finite time. In their 
investigations and other studies on disentanglement in 
open quantum systems [6] [7], only qubit systems are 
considered. Here, the disentanglement process is char- 
acterized by time evolution of the concurrence [8]. It is 
well-known that the concurrence is not available in the 
higher-dimensional systems. 
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For the systems with spins larger than 1/2, one can 
use the positive partial transpose (PPT) method [9] to 
study disentanglement. For the mixed states of two spin 
halves and (1/2 ,1) mixed spins, the PPT method can 
fully characterize entanglement. However, in the case of 
two qutrits and even larger spins, one only know that if 
a state does not have PPT, the state must be entangled. 
In other words, one may use the method to witness en- 
tanglement. Actually, in the case of higher dimension, 
there are two qualitatively different types of entangle- 
ment [10], free entanglement (FE) which corresponds to 
the states without a PPT, and bound entanglement (BE) 
corresponding to the entangled states, however, with a 
PPT. The BE is an intrinsic property and cannot be dis- 
tilled to a singlet form, thus it cannot be used alone for 
quantum communication. Nevertheless, the BE can be 
activated and then contribute to quantum communica- 
tion [11]. A formal entanglement-energy analogy [12] im- 
plies that the bound entanglement is like the energy of 
a system confined in a shallow potential well. If we add 
a small amount of extra energy, behaving as a perturba- 
tion, to the system, its energy can be deliberated. The 
existence of bound entangled states reveals a transparent 
form of irreversibility in entanglement processing [13]. 

In this paper, we consider an open composite system, 
a two-qutrit system commonly coupled to an environ- 
ment, and study a type of dynamical process of disen- 
tanglement, where the two qutrits are initially prepared 
in a bound entangled state. We would like to reveal 
that different environment gives different dynamics of en- 
tanglement. Firstly, a bosonic heat bath is considered. 
We remark that this modeling of environment is univer- 
sal [14, 15] in the sense that any environment weakly cou- 
pled to a system can be approximated by a collection of 
harmonic oscillators. Secondly, we consider a spin envi- 
ronment consisting of spin halves which can be considered 
as a fermionic environment. We let two types of bound 
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entangled states being initial state of the two qutrits in 
order to find the different properties of bound entangled 
state during the quantum dephasing. And initially the 
environments are assumed to be at thermal equilibrium 
states, which helps us to find effects of the thermal fluc- 
tuation on dynamics of quantum entanglement. 

This paper is organized as follows. In Sec. II, we con- 
sider the bosonic environment and give the analytical 
results of FE and BE. we numerically study the BE to 
illustrate the details of the dynamics of entanglement. 
In Sec. Ill, the two qutrits are coupled to a spin environ- 
ment. Also the analytical and numerical results are given 
to show the effects of coupling strength, temperature and 
energy spectrum structure on the dynamical behaviors of 
entanglement. The conclusion is given in Sec. IV. 



II. BOUND ENTANGLEMENT IN A BOSONIC 
ENVIRONMENT 

We start with a well known model of the pure dephas- 
ing [15, 16], where two qutrits interact with the envi- 
ronment, which is modelled as a heat bath with many 
harmonic oscillators of frequency u>j. The model Hamil- 
tonian reads 



H = ^2 H 3= fl^ b h + 9(b] + bj)(Sn + S 2z )], (1) 



is the displace- 



and D(z S z ) — exp (zfi^ — z*b^j Sz 
ment operator. 

Before discussing the dynamical process of entangle- 
ment, we introduce two quantities to quantitatively study 
entanglement. One is the negativity [18], which can be 
used to study FE. For a state p, negativity is defined in 
terms of the trace norm of the partial transposed matrix 



N( fi ) = 



Ip Ti IIi-1 



(7) 



where T\ denotes the partial transpose with respect to 
the first subsystem. If M > 0, then the two-spin state 
is free entangled. As an entanglement measure, the neg- 
ativity is operational and easy to compute, and it has 
been used to characterize entanglement in large spin sys- 
tem very well [19]- [21]. 

In order to characterize BE, one can use the so-called 
realignment criterion (cross-norm criterion) which proved 
to be very efficient [22]. The operation of realignment on 
the density matrix is just as (p R )ij.ki — Pikji- A separate 
state p always satisfies \\p R \\ < 1. Thus, a quantity for 
the BE can be defined as 



K(p) = max{0, \\p R \ \ - l} . 



(8) 



We call this the witness quantity. Only when IZ(p) > 
and N{p) =0, the state is bound entangled. 



where &j and bj are creation and annihilation operators, 
respectively, S\ z and S 2z are z components of two spin-1 
operators, and g denotes the coupling strength between 
the spins and the heat bath. 

In order to study the dynamical process of entangle- 
ment in our system, it is convenient for us to study the 
time evolution in the interaction picture. Here, 



H = ^ hujjbjbj 



(2) 



A. Horodecki's Bound entangled state 

In the following discussions, we consider the dynamical 
evolution process of the two-spin 1 system, deriven by 
the Hamiltonian (1) with the initial state being in the 
Horodecki's bound entangled state [11]. 



1. Analytical results 



is the free Hamiltonian, and the interaction Hamiltonian 
Hi = Y,a(b] + b ] )(S lz + S 2z ). (3) 



The bound entangled state reads [11]: 

< . 2 n a 5 — a 
P„(0) = -P+ + -Q+ + —6- 

2 < a < 5, 



(9) 



Then, through the Wei-Norman method [17], the time 
evolution operator in the interaction picture is factorized 
as, 



U(t) = Y[e i *^ s *D[z j (t)S z } 



where S z = S\ z + S 2z , 
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(1 - e*"'*) 



(4) 

(5) 
(6) 



where 



p+ = |* + )<*+|,|* + ) = 

-L(|00) + |ll) + |22)), 
Q+ = i(|01)(01| + |12)(12| + |20)(20|), 



(10) 



Q- = g(|10)(10| + |21)(21| + |02)(02|). (11) 

where \m1m2) , (mi, m 2 = 0,1,2) are the eigenvectors 
of S z = Si z + S 2z , with the corresponding eigenvalues 
mi + m 2 — 2, respectively. 
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In Rcf. [10], Horodecki demonstrated that 

separable for 2 < a < 3, 
p a is ^ bound entangled for 3 < a < 4, 
free entangled for 4 < a < 5. 



(12) 



And the density matrix for the initial state of the total 
system is a simple direct product 



Ptot (°) =Pa®PEi 



(13) 



where p E is the density matrix of environment. 

Driven by the time evolution operator (4), the system 
will evolve from the bound entangled state p a into the 
state described by 



Pl,2 (*) 



Tr E [U(t)p tot (0)f/t (t) ] 

2 r (|oo>(oo| 



11) (11| + |22) (22|) 



where 



21 

+ (Fi(t)|00)<ll|+H.c.) 
+ (F 2 (t)|22)(ll| + H.c.) 
+ {F 3 (t)\00) (22| + H.c)] 
a 5 — a 



Fx (t) = Tr s [ PE U\ (t) U (i) 
F 2 (t) = Tr s \p E U\ (t) U 2 (t) 
F 3 (t) - Tr E \p E u\ (t) U (t) 



(14) 



(15) 



are decoherence factors [15] . The unitary operators Uq (t), 
Ui (t), and XJ<i it) are derived from Eq. (4) just by replac- 
ing operator S z = S\ z + S2 Z with numbers —2,0 and 2, 
respectively. 

From the reduced density matrix (14), the realigned 
matrix becomes 



A,^,^ © © -^2x2 © -^2x2 



)• 




1,2,3). 



Then, the witness quantity 1Z is obtained as 
U(p h2 ) = max{||pf 2 (t)||-l,0} 



(16) 



= — max{\/3a 2 - 15a + 19 



+2(|*i| 



\F, 



\F* 



7,0}. (17) 



As mentioned above, the positive witness quantity can 
quantify the nontrivial BE only when the negativity van- 
ishes. Thus, we need to calculate the time evolution of 
negativity. 



We first make the partial transpose of p 12 with respect 
to the second system and obtain 

{Pl2{t)) T2 = TTf (^3x3 © D 2x2 © ^2x2 © ^2x2) , 




D. 



(fe) 

2x2 



(2P* 5-a) ( fc - !' 2 ' 3 )- 



(18) 



Then, from the above equation, we immediately obtain 
the negativity 

1 3 I 

jv(pi, 2 ) = 42 E max {°- v ( 2a - 5 ) 2 + 16 i^i 2 - 5 1- 

fe=l 

(19) 

Thus, we have obtained analytical expressions of witness 
quantity 1Z and negativity Af in terms of the three de- 
coherence factors. It is natural to see that if the de- 
coherence factors are zero, namely, the completely de- 
coherence occurs, from Eqs. (17) and (19), we have 
1Z = Af =0. From Eq. (19), we can also see that in the 
region 3 < a < 4, negativity always gives zero at any 
time since \Fk\ < 1. 

From the above discussions, once we know the deco- 
herence factors, the two quantities 1Z and Af for detect- 
ing entanglement can be determined. So, we are left to 
obtain these decoherence factors. It is well known that 
high temperature may enhance the decoherence, thus it 
is reasonable to choose a thermal equilibrium state as the 
initial state of the heat bath, which is described by the 
density matrix 



Pe = 



Tr e 



-f3Hujjbj bj 



= J|(l - e ~ phui ) e -P h "i b ] b i j 



(20) 



where (3 — 1/fcgT, k E is the Boltzmann's constant, and 
we choose k E — 1 for simplicity in the following. 

For the bosonic environment we calculate the decoher- 
ence factors in the coherent-state representation. The 
P-representation for the thermal state is given by 

Pe (0) = Hpej=Y[J Pj («) |«> ("I d 2 a, (21) 

(22) 



1 I \a 

Pi \ a ) = _ - ; ex P 



Tr (nj) 



( n j) / ' 



-l . 



where (rij) = ( y e f3huJj — l) "is the thermal excitation 
number of harmonic oscillators. From Eq. (15) and us- 
ing the P-represcntation, one obtains the modulus of the 
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FIG. 1: The modulus of the decoherence factor |-Fi(t)| versus 
time with g — 2. The frequencies of heat bath are chosen 
randomly in a region ujj £ [50,55]. The size of the bath is 
L — 200 and the system is at the temperature T = 1. The 
horizontal line in the figure corresponds to the threshold value 
\Fi\ = 0.8398. 



decoherence factors decay as a Gaussian or a exponen- 
tial form with time, and consequently we know the time 
behaviors of the BE. 

i) Let us choose a certain type of distribution of u>j 
in the region [0, u] , where oj is an arbitrary value larger 
than zero. We do not care about the exact form of the 
distribution, however it can be achieved for us to choose 
a sufficiently small cutoff frequency ujj c to make sure that 
at a finite time, the decoherence factor (h = 1) 



\F x (t)\ < exp£ 



Z^sm 2 f^)(2(n,) + l) 



exp 



-2 5 2 ^(2(n J ) + l)t 2 



-7f 



(27) 



decoherence factors [231 



where 



|*i (*) I = 1*2 (*)| = J] l Tr ^ [p Ej D(2z j )}e' 

3 

= TJ e -2| 2:3 | 2 (2<r l3 > + l) 

j 

sg 2 



1*3 W 



4<I> j 



ex P(ElS(2(^) + l)sin 2 (^)i2:; 1 



Y[Tr j (p Bj D(-4z J ))=\F 1 (t) 



(24) 



1) 



(28) 



It can be seen that the decoherence factor displays a 
Gaussian decay with time. Moreover one may observe 
that the decay parameter 7 increases at high tempera- 
ture since (nj) is a monotonically increasing function of 
temperature T, and enlarging the strength g can also 
increase 7. Substituting Eq. (27) to (25) leads to 



As expected, the above three quantities are smaller 
than or equal to unity. Now, we study the decoherence 
of BE, and choose parameter a = 4 in the bound en- 
tangled state in the following discussions. This choice 
of parameter maximize the quantity 1Z. Then, Eq. (17) 
simplifies to 

npia) = ^ max {o, 2 (2 \F X \ + |^i| 4 ) + V7 - 7} . 

Then, we find that the dynamic properties of BE is 
thus directly related to the one single decoherence fac- 
tors |*i(i)|. By numerical calculation, one obtains the 
threshold point of 



LFjl « 0.839829, 



(26) 



before which the quantity 1Z is larger than zero, im- 
plying that the state is a bound entangled state. In 
Fig. 1, we numerically show the modulus |-F\(i)| versus 
time. The frequencies uij are chosen randomly in a region 
ujj € [50,55]. Then, the modulus |-Fi(i)| oscillates with 
time and periodically crosses the horizontal line corre- 
sponding to the threshold value Fi = 0.8398. Obviously, 
the witness quantity displays discontinuous behavior and 
below the line it becomes zero. 

In the following we consider some special cases of the 
energy distribution in the environment and find that the 



ft(Pi, 2 ) = 7T7 max {o, 2 + e^) + Vl - l) , 

and it will decay to zero in a fixed time to, which can be 
determined from Eq. (26) 



t = 0.4176/V7- 



(30) 



When the evolution time is larger than the threshold 
value to, the BE suddenly vanishes. 

ii) If we choose some continuous spectrum, the sum in 
the decoherence factors (we assume Ti = 1) becomes 



ln|*i(t)| = -X; 



85, . 2 / UJjt 



to 



3 



(31) 



where gj - 
tion p (u>j 



1). Assume a spectrum distribu- 



In |*i (t) I 



, the above equation becomes 



■ 2 , 
sin —z—du)4. 
2 J 



(32) 



For some concrete spectrum distributions, interesting cir- 
cumstances may arise. For instance, when p(u>j) = 
r y/(2TrgJ) the integral converges to a negative number 
proportional to time t, precisely, \Fi (t) | = exp(— jt), 
1*3 (t) I = exp(— 4r/t). Thus, in this case, the reasonable 
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FIG. 2: TZ versus time with different coupling parameter g. 
The frequencies of heat bath are chosen randomly in a low 
region u)j £ [0, 5]. The size of bath L — 200 and the system 
is at a finite temperature T = 1. 




FIG. 3: (a) TZ versus time with different coupling parameter 
g. The frequencies of heat bath are chosen randomly in a 
higher region u)j £ [50, 55]. The size of the bath L = 200 and 
the the 
freq iling 
9 = 



assumption on the energy distribution brings us a ex- 
ponential decay of decoherence factor and entanglement 
with time. 

m) Now we will choose another more general distri- 
bution p(u>). Assume that all the coefficients gj are 
equal: gj = G. If the frequencies lie within an interval 
[a; 2] and the distribution is homogeneous, we have 
p («*) = N/ [u 2 - wi) , thus [24] 



ln|Fi(<) 



E 

3 



8g 2 



f sin 2 ^1 
2 9 



sin 2 (cjfci) 



ckvi 



< 



-2G 2 N 
uj 2 — 

-2G 2 N 4 



2 sin u>kt 



uj 2 - 
—2G 2 N 



duj k 

2 

sin 2 (cjkt) duik 

2 

2 cos 



^±^t) sin I 



(uj 2 - LOl) t 



.(33) 



By substituting the above equation into the Eq. (25), we 
see that when the environment has sufficiently large size 
L, the decoherence factor and the quantity TZ will decay 
with time rapidly. 



2. Numerical results 

Next, we resort to numerical calculation to test the 
above analysis and show more dynamic behaviors of en- 
tanglement. 

In Fig. 2, we choose a random distribution of the envi- 
ronment energy Wj over a finite region [0, cj], and give the 
time behaviors of the quantity TZ. A Gaussian decay is 
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\ T=50 




\ T=100 








FIG. 4: (a) Three dimensional (3D) diagram of TZ versus 
time and temperature, with L — 200 and g = 3. The region 
of frequencies is uij € [50, 58]. (b) Quantity TZ versus time at 
different temperatures, (c) TZ versus temperature at a fixed 
time t = 0.003 for different sizes. 



exhibited, and this is consistent with the previous anal- 
ysis (Eq. (27)). Of course, larger g accelerates the decay 
process of the BE. When we change other parameters, 
such as the width of frequencies and the temperature, 
the quantity TZ displays similar behaviors as long as the 
region of ujj begins from zero. This implies that the har- 
monic oscillators with lower energies in the heat bath 
determine the behaviors of the BE. . 

Instead of choosing the low energies of the harmonic 
oscillators in environment, we consider a random distri- 
bution of ujj in a higher frequency region [50,50 + 5], 
where S is the width of the distribution. The numerical 
results are shown in Fig. 3. For small values of the cou- 
pling constants, as shown in Fig. 3 (a), the BE displays 
oscillations with time. For g = 2, we observe collapses 
and revivals of the BE. The revivals result from the re- 
vivals of the decoherence factors. When the coupling 
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strength is strong enough, the BE decays rapidly to zero 
without revivals. 

Fig. 3(b) presents numerical results for different fre- 
quency widths. Consider an extreme case 5 = 0, that 
is, only one frequency is taken into account. In this case, 
the BE displays periodic collapse and revivals for large g, 
and the revival amplitude of 1Z is one. From the figure, 
we see that when the frequency width increases, the re- 
vival amplitude decreases. Increasing the width S means 
that the harmonic oscillators in the heat bath own much 
more different frequencies and the entanglement revival 
will be suppressed. 

Now we consider the thermal effect on the BE in Fig. 4. 
The subfigures (a) and (b) show that the thermal fluctu- 
ation can destroy entanglement and accelerate the decay- 
ing process. With the joint effect of thermal fluctuation 
and strong coupling g, entanglement vanishes in a finite 
time without reviving. In Fig. 4(c), we can see that with 
the temperature increasing the entanglement decreases 
to zero, and enlarging the size of heat bath can suppress 
the BE and accelerate the decay. 



B. Second bound entangled state 

We choose another 3x3 bound entangled state as the 
initial state of the two qutrits which was introduced by 
Bennett et al. [25] from the unextendible product bases: 
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FIG. 5: We consider the BE and FE together. And in the 
four subfigures the blue lines with circle markers correspond 
to 1Z and the red lines with x markers denote jV. (a) 1Z and 
j\T versus time with coupling g = 1, the size of environment 
L = 300, and the system is at temperature T = 10. The 
frequencies of heat bath are chosen randomly in the region 
ujj G [50, 55]. (b) Only changing the coupling to a larger one 
g = 5, and the other parameters are the same as subfigure 
(a). In the figure (c), parameter 8 = 9, the other parameters 
are the same as subfigure (b). In subfigure (d), it shows 1Z 
versus temperature at t = 0.005 and j\f versus temperature 
at t = 0.115. g = 1 and L = 300. 



4l°)d°) - I 1 )). 
i(|0)-|l»|2>, 

^|2>(|1> - |2», 

i(|l))-|2))|0), 

i(|0) + |l) + |2»(|0) 



(34) 



ID 



|2», 



from which the density matrix could be expressed as 

4 



i 4 

9 = 1(^x9-^1^X^1), 



(35) 



3=0 



In this case, the dynamics of entanglement are deter- 
mined by six decoherence factors, and analytical results 
are not available. We numerically calculate BE and FE, 
and the results are shown in Fig. 5. 

In order to compare BE and FE, we have numeri- 
cally given the time behaviors of both quantity 1Z and 
negativity M . We choose a higher frequency region 
ujj G [50, 50 + S], which will induce some interesting prop- 
erties of BE and FE. In Fig. 5 (a), we see that the nega- 
tivity can be nonzero, in contrast with the first BE given 
by Horodecki et al. The nonzero negativity implies that 



the state is free entangled. The state keeps bound entan- 
gled for a short time, and then, the state is free entan- 
gled. When the coupling g becomes stronger, as shown 
in Fig. 5 (b), after some revivals, both quantities become 
zero. 

In Fig. 5 (c), we extend the width of frequency region 
to 5 = 9. Expanding the energy region makes the par- 
ticles in heat bath own more chances to take different 
energy. And this will prevent the revival of entangle- 
ment. In Fig. 5(d), we show the behaviors of 1Z and j\f 
against temperature for a fixed time. With the increase 
of temperature, quantity TZ and negativity Af decrease 
gradually, and at last the thermal fluctuation destroys 
the entanglement completely. 



III. BOUND ENTANGLEMENT IN A SPIN 
ENVIRONMENT 

To show the universality of the influence of environ- 
ment on the time evolution of BE, we need to use differ- 
ence modeling of environment. Here, we consider an envi- 
ronment consisting of N spin halves. The corresponding 
model Hamiltonian reads [26] 



H 



fJ 



(Si 



fe=l 



•(*) 



(36) 



7 



(k) 

where a \ denotes the z-componcnt of the Pauli vec- 
tor, and g denotes the coupling strength between central 
spins and environment. We notice that the above model 
has been considered by Zurek [26] as a solvable model of 
decoherence. 

The time evolution operator can be expressed as: 



U(t) = l[exp(-itAu k a^), 



(37) 



fc=i 



where we define a special operator-valued parameter 



A = | (S lz + S 2z ) . 



A. Horodecki's Bound entangled state 



(38) 



In a similar vein as the discussions of entanglement in 
the bosonic environment, we first study the disentangle- 
ment of Horodecki's bound entangled state and give the 
analytical results. 



1. Analytical results 

Let us consider the whole system initially starts from 
a product state 

Plot (0) =Pa®PE, 

where the initial state of the two qutrits p a is a mixed 
BE state represented in Eq. (9), and p E denotes the ini- 
tial state of the environment which is assumed to be a 
thermal state described by the density matrix 



p E = U 



CO 



L 2cosh(/3cjfc) 



(39) 



Then the reduce density matrix at time t can be given 
by the same matrix as Eq. (14), and now the three de- 
cohrence factors in this spin environment can be obtained 
as 



fc=i 



\Fs(t)\ = 




k=l 



sin 2 (gtuj k ) 
cosh 2 (/3wfe) 



sin (2gttOk) 
cosh 2 (/3wfc) 



(40) 



(41) 



From the Eqs. (40) and (41) one can find each deco- 
herence factor | Fk | is less than unity, which implies that 



in the large L limit, \Fp. (t)\ will go to zero under some 
reasonable condition. Now, we make some further anal- 
ysis by introducing a cutoff number K c similar to the 
discussion in Ref. [27]. We define the partial product as 



l^i (t)\ c 



U\Fi,k\ 

fc>0 



>l*i(*)l 



(42) 



from which the corresponding partial sum ln|fi (t)\ c = 

~ Sfc>o l m ^i,fc|- We can do some heuristic analysis in 
some special conditions such as confining the energy spec- 
trum in a region from zero to a nonzero value Uk <G [0, u]. 
When the cutoff number K c is small enough, in a finite 
long time, with some proper g we can pick out some 
tiny u)k to make gtuk begin a small one and achieve the 
approximation sin 2 (gtujk) ~ (gtwk) 2 - At a finite temper- 
ature we can have 



H*i(*)| c = oE ln 1 



(gtu> k ) \ 



cosh 2 (/3wfe) / 



fc>0 

iff, uj ' 
2 I ^ cosh 2 (0u>k) 



9 2 t 2 



\fc>o 



- 7 t 2 



(43) 



where 



7 = 77.9' 



E 



cosh (/3wfc) 



From Eqs. (42) and (43), we find that the decoherence 
factors decay in a Gaussian form with time, therefore 
from Eq. (25) it is apparent that the witness quantity 
1Z will vanish in a finite time. Also from (43), if the 
temperature is very low and quite nearly to T = 0, 7 
approaches zero, and the decoherence factors quantity TZ 
will not decay with time. It implies that, in our system, 
the temperature greatly affect the dynamics of BE. It is 
a rough calculation in our analysis, nevertheless it gives 
us some constructive results. 



2. Numerical results 

If the frequency distribution are in the region LUk G 
[0, u], the decoherence factors displays a Gaussian decay, 
which is analytically studied in the former section. And 
in Ref. [28], Gaussian decay was shown numerically in 
various distributions of couplings. So, here, we consider 
the distributions such as [50, 50 + 5], a higher frequency 
region. In Fig. 6 (a), the larger coupling strength g makes 
a stronger oscillations of BE, which is different from the 
case of bosonic environment. Mathematically, we can 
understand that from Eq. (40) and (41), parameter g 
can change the frequency of the periodic function. The 
collapses and revivals of 1Z is also observed here, and 
the revival amplitude decreases with time. Fig. 6 (b) is 
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FIG. 6: (a) TZ versus time with different coupling parameter 
g. The couplings uik of the spins in environment is random 
in l nt 
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FIG. 7: (a) Three dimensional (3D) diagram of the TZ ver- 
sus time and temperature, with L = 300, g = 0.5 and 
uJk £ [50,55]. (b) Several part sections of the 3D diagram. 
Explicitly, it presents TZ versus time at different temperatures 
of T = 1, T = 10 and T = 40. (c) shows TZ decays with 
temperature at a fixed time t = 0.005 with different sizes of 
environment L = 300, 1 000, and 5000. 
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FIG. 8: TZ versus time at different temperatures, T — 10, 
T = 15 and T = 35. With the coupling g — 0.5, environment 
size L = 300 and lu^ £ [50, 55]. 



stable against time for low temperatures. This is a quite 
different property from the case of bosonic environment. 
In Fig. 7 (c), we also plot the BE against temperature 
for different size of environment. As we expected, the 
larger size of the spin environment accelerates the decay- 
ing process. 



B. Second bound entangled state 

We now consider the case that the two qutrits initially 
in the second bound entangled state (35) in the spin en- 
vironment. We choose cut G [50, 55] and study the finite- 
temperature effects on entanglement. In Fig. 8, we can 
see the similar phenomena with Fig. 6, namely, at low 
temperature BE only oscillates around the initial value 
of the BE. At higher temperature BE decays sharply with 
time, and revivals are also observed with revival ampli- 
tudes decreasing with time. We also study negativity, 
however, it always keeps zero which means there does 
not exist FE all the time. 



a plot of TZ for different width of frequency distribution. 
The wider distribution will smear the collapse and revival 
phenomenon, namely, the BE evolve smoothly with time, 
and the BE is always there. 

Effects of the thermal fluctuation on the dynamic of 
BE are shown in Fig. 7. The 3D plot (a) displays a 
flat at low temperatures about T < 10, which implies 
that the BE is stable against temperature in this region. 
When temperature is high enough, BE rapidly decreases 
to zero. In Fig. 7(b), we can see explicitly that temper- 
ature plays an important role in the dynamics of BE in 
this spin environment. It is just like a control process that 
when temperature is higher than some value, BE will de- 
cay sharply with time. At a very low temperature, from 
Eq. (40) and Eq. (41), we know that decoherence factors 
are approximately one with very small oscillations. And 
the BE is of a little change with time, namely, the BE is 



IV. CONCLUSION 

In summary, we have studied the decoherence phenom- 
ena of two-qutrit system couple to an enviroment when 
this open systems are initially prepared in a bound entan- 
gled state. The so-called " sudden death" phenomenon [5] 
of entanglement can be found in our studies as a com- 
mon feature of BE evolution. Two typical pure dephasing 
systems, the bosonic and the spin systems are considered 
in order to show this kind of universality of decoherence 
process with BE. Here, we used the realignment crite- 
rion to characterize BE, and the PPT criterion to study 
FE. Beyond the negativity, we have introduced a novel 
witness quantity TZ to study BE. Those two approaches 
are operational and convenient to use. Two kinds of BE 
in initial state of the open system are considered, one is 
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given by Horodecki et al, and another is constructed from 
the unextendible product basis. 

One of our central result is to express the quantity 
1Z and negativity Af in terms of three decoherence fac- 
tors, and these factors are analytically obtained. In the 
case of bosonic environments, the Gaussian decay and 
the exponential decay of the BE was found, and in the 
case of spin environment, we find that BE can display a 
Gaussian decay. In both environments, the collapse and 
revivals of the BE are observed for system frequency dis- 
tribution being in a higher region with an appropriate 
width. Larger coupling strength g, larger environment 
size L and higher temperature will enhance the disen- 
tanglement process. For the Horocecki's BE in the spin 
environment, we find that the BE can be stable against 
low temperature increase. 

Finally we have to point out that since we are lack of 
entanglement measure for two qutrits, the study of de- 
coherence of entanglement here is incomplete. Neverthe- 



less, the realignment criterion and the PPT criterion are 
very efficient to characterize BE. It will be interesting to 
consider decoherence of BE under other decoherence pro- 
cesses such as dissipation, and investigate the robustness 
of the BE. 
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